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Abstract 

We report a method for probing electromechanical properties of multiwalled carbon nan- 
otubes(CNTs). This method is based on AFM measurements on a doubly clamped suspended 
CNT electrostatically deflected by a gate electrode. We measure the maximum deflection as a 
function of the applied gate voltage. Data from different CNTs scale into an universal curve within 

o, 

the experimental accuracy, in agreement with a continuum model prediction. This method and 
the general validity of the scaling law constitute a very useful tool for designing actuators and in 
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£Sj ■ general conducting nanowire-based NEMS. 
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Carbon nanotubes (CNTs) are promising candidates for designing and developing nano- 
electromechanical systems (NEMS) because they combine excellent electronic and mechani- 
cal properties. High conductivity of CNTs allows for designing simple sensing and actuation 
systems based on the direct electrostatic coupling with metallic gates. Their exceptional 
stiffness, low mass and dimensions ensure operating frequencies in the GHz range making 
them suitable for a number of applications^ |. Some prototypes of CNT-based NEMS such 
as memory devices Q], nanotweezersj^, high frequency oscillators j^J and actuators have 
already been demonstrated. Theoretical studies of CNT-based switches have been recently 
published [30. Particularly interesting are molecular dynamics simulations of Dequesnes 
and coworkers (J which have made evident the validity of continuum models (beam theory) in 
describing CNT deformation when its length/diameter ratio is larger than 10. This is a very 
important point because it simplifies the description of CNT mechanical deflection avoid- 
ing expensive and time consuming atomistic simulations. However, designing CNT-based 
NEMS requires a very precise knowledge of the static deflection under actuation. Since in 
most of the CNT-based NEMS actuation is electrostatic the relevant parameters are: the 
geometry of the CNT (inner and outer diameter, length), its Young's modulus Y, and the 
geometry of the device which conditions the CNT- gate(s) electrostatic coupling and thus 
the actuation efficiency. The interplay between the different geometrical parameters indeed 
makes theoretical prediction tools indispensable to properly scale any practical device based 
on suspended CNTs. Comparison with experiments is equally important both to evaluate 
the validity of the predictions and precisely measure Y, since it determines most of the static 
and dynamic properties of CNTs. 

In this letter, we report an on-chip test method for measuring the deflection of suspended 
and electrostatically actuated CNTs. We develop a theoretical framework for the modeli- 
sation of the electromechanical behavior of CNTs based on continuum models and validate 
it experimentally u sing that test method. Furthermore, we determine the CNT- Young's 
modulus [^J 1]] [l4|0] very precisely. The heart of the system is a doubly clamped 
suspended CNT three-terminal device (see Fig. 1) deflected by an electrostatic force induced 
by a back gate. In this geometry, Van der Waals forces (| can be neglected. The CNT de- 
flection only depends on its physical properties and on the electrostatic environment given 
by the connecting and gate electrodes, thus allowing a quantitative comparison with calcu- 
lations. We measured the maximum deflection Umax (umax = H — y(L/2) ) see Fig. la) of 
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a suspended multiwalled CNT as a function of the applied gate voltage Vq. We show that 
data from different CNTs (with different diameters D and lengths L ) scale into an universal 
curve within experimental uncertainties. This scaling law that we derived from continuum 
beam theory when Umax *C D, can be extended to the Umax > D range, where the induced 
stress T due to elongation of the CNT becomes important. Furthermore it is also valid in 
the general case where there is an electrostatic force profile applied along the CNT. The set 
of consistent data allows us to accurately determine Y = 0.41±0.05 TPa. Our method and 
the general validity of the scaling law constitute a very important outcome for designing 
conducting nanowire-based NEMS. 

A schematic diagram and a typical SEM picture of a CNT-based device is presented in 
Fig. 1. A multiwalled CNT of diameter D is clamped by two metallic pads and suspended 
over a length L on top of a highly doped silicon substrate that acts as a gate. The distance 
between the CNT and the gate is fixed by the sacrificial silicon dioxide layer thickness 
H (230 nm in the present case). The multiwalled CNTs used were synthesized by arc- 
discharge evaporation and carefully purified to remove amorphous carbon and graphitic 
nanoparticles 14( . The sample were fabricated as follows: some droplets of a sonicated 
suspension of CNTs in dichloroethane (50 fi g/ ml) are deposited on an oxidized Si wafer, 
and blown dry after one minute under nitrogen flow. Atomic force microscopy[15] is used 
to image, select, and locate CNTs with respect to pre-patterned alignment marks. Selected 
CNTs are then connected with two electrodes designed by electron beam lithography, and 
deposited by thermal evaporation of gold (70 nm) with a chromium adhesion layer (0.5 nm). 
These metallic electrodes are used both as conducting electrodes and anchor pads. The 
sample is dipped in BHF to remove all the Si02 underneath CNTs, rinsed in DI water and 
ethanol and dried on a hot plate at 50° C. 

The sample is mounted in a commercial AFM equipped with a conducting tip. CNTs and 
AFM tip are grounded, in order to avoid any electrostatic interaction between them [la]. To 
measure the deflection of the CNT when a voltage Vq is applied to the back-gate electrode, 
the AFM is used in the tapping mode. The ex per imental setup is depicted in the insert of 
Fig. 2. The AFM tip is placed and immobilized [l7|] at the center of the CNT. The feedback 
voltage applied to the piezoelectric element V P i ezo , that controls the vertical position of the 
tip, is monitored as a function of Vq. Fig. 2 shows a typical result for a suspended CNT, 
600 nm long and 10 nm in diameter (device #1). As \Vq\ increases, the CNT deflects 
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downwards, the AFM tip has more room to oscillate and the oscillation amplitude increases. 
To keep constant the oscillation amplitude of the tip, a positive V P i ezo is applied by the 
AFM electronics. The measured V P i ezo , proportional to the deflection of the CNT, is found 
to vary as Vj. This is expected for a CNT that follows the Hooke's law, because the applied 
force is proportional to Vj (see EqJ3] below). The AFM tip position is also influenced by 
the electrostatic attraction to the substrate. Indeed, when the tip is placed on the sacrificial 
silicon dioxide layer at the same CNT-gate distance (position 2 in Fig. 2), increasing \Vq\ 
increases the attractive electrostatic force between the tip and the substrate, the oscillation 
amplitude decreases and a negative V P i ezo is applied by the feedback loop. The genuine 
CNT deflection is obtained by subtracting from the curve measured at position 1 the one 
measured at position 2 previously divided by the relative permittivity of the silicon dioxide 
(dotted line in Fig.2)|]j|. 

We repeated the experiment on four devices with different L and D parameters [19]. Using 
the calibration of the piezoelectric element j^J (243.8 nm/ V) the deflection of the CNT, 
Umax, can be obtained as a function of the applied Vq- Fig. 3a shows results for three 
different devices: #1 (L = 600 nm, D = 10 nm ), j^l {L = 770 nm, D = 13.5 nm ) and 
#3 (L = 480 nm, D = 10 nm ). The parabolic behavior is made evident in this Log-Log 
scale representation (the dotted line is a guide to the eye that has Vq dependence). 

The strong dependence of Umax (Vq) on different CNT parameters (L and D) can be 
understood using a continuum beam equation |2]J for describing the deflection of the CNT 
(denoted by y(x)) when an electrostatic force per unit length F e i ec is applied. In reduced 
units (y — y/D and x = x/L) this equation reads 

A TL 2 d 2 y L 4 

dx~* ~ YTdx^ = yJd ec ^ Vg) (1) 

Here Y is the Young's modulus (ty pic ally of the order of 1 T Paj9j-[l3|) and / the moment 
of inertia of the CNT approximated |22| by / = 7rD 4 /64. T is the stress force induced by 
elongation of the CNT when the electrostatic force F e i ec (y, D, Vq) is applied. T can be 
calculated by, 

nD 4 Y f 1 (dy x 2 



T = — 77T- l-r) dx (2) 
8L 2 Jo \dxj K ' 

To solve Eq^, we have to assume that the stress force T has a constant value, which is ob- 
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tained later from the self-consistent condition (Eq|2J). The electrostatic force F e i ec (y, D, Vq) 
caa be m ode,ed by ordering the electrostatic ene rgy E aec = JC D) V 3 and the prin- 
ciple of virtual work|23J], C (y,D) being the capacitance per unit length of a CNT over a 
metallic plane. This force has the form: 

F elec (y,D,V G ) = ^0g(2y) (3) 



Where Eq is the permittivity of vacuum and g (t) = Jt (t + 2) In ( 1 + t + Jt (t + 2) 



This function can be well approximated in the range of practical interest |24| by g (t) = At^ 13 
with A = 0.270 and (3 = 1.456 . Since in our case H ^> Umax, considering the electrostatic 
force (at a given Vq) as a constant is a good approximation. This is strictly true for an infinite 
geometry. Indeed, contacting electrodes modify the electrostatic force near the CNT-metal 
interface due to screening and the AFM tip (in spite of being at the same potential as the 
CNT) modifies the electrostatic force near its position. We calculated the electrostatic force 
profile that builds up by 3D numerical simulations (FEMLAB) of our structures. It turns 
out that the deflection of the CNT produced by the electrostatic force profile is equivalent to 
the one generated by a constant force F e i ec (H, D, Vq) times a correcting factor. This factor 
C FP , is: 1.26 for device #T 1.28 for #2, 1.45 for device #3 and 1.21 for device #4 ( L 
= 730 nm, D = 15 nm) [2^j. 

It is interesting to note that the second term in EqQis negligible}^ at low Vq, where 
Umax <C -D. In this limit the maximum deflection is well approximated by 

Umax T ^fV G L 2 \ 2 jz C FP Ae (2Hy f3 

— = K {vrm) > K = w (4) 

Fig. 3 depicts the results on the four devices rescaled using Eq|U The curves indeed 
collapse onto a single curve within the experimental accuracy in determining L and D. 
This finding implies that Y is nearly the same for all the devices investigated, in agreement 



with theoretical predictions |27l|. Furthermore, it allows a quantitative determination of 
the Young's modulus as shown in the insert of Fig. 3b, where the points represents Y 
calculated from the data using EqJ3]for the four devices investigated. The average Y value 
obtained from these experiments is 0.41 ± 0.05 T Pa. This value is smaller than expected 
compared to previous measurements, although it is statistically indistinguishable from the 



one obtain by J. P. Salvetat et al. |12j on the same source of CNTs. One possible explanation 
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might be defects induced either by the purification treatment or by the lithographic steps. 
Experiments are underway to examine these issues. 

In order to explore the validity of the scaling law predicted and experimentally observed at 
low deflections, we solved Eq^on a larger Vq range. Fig. 4a shows the calculated maximum 
deflection of the CNT [umax) as a function of Vq for the structure depicted in Fig. la. 
Different curves correspond to different D, L parameters (experimentally accessible) for a 
fixed H value. Since H is still much greater than Umax, considering the electrostatic force 
along the CNT as a constant remains a good approximation. Notice the large variation of 
Umax (Vq) when L or D varies. This simply reflects the strong dependence on geometrical 
parameters. Now, if the curves depicted in Fig. 4a are rescaled according to Eq0] a collapse 
onto a single curve is made evident (see Fig. 4b). This scaling law, expected at low Vq 
as shown above, proves to be valid in the whole range investigated. This striking result 
observed for a constant force per unit length can be generalized to the case of multiple gates 
as long as the functional form of the force profile along the CNT is the same for devices with 
different L and L>|3]- 

In conclusion, we have presented an accurate method for probing the electromechanical 
properties of suspended nanotubes and measuring the Young's modulus using a conducting 
AFM on a simple test bench structure. In this method the maximum deflection Umax of the 
CNT is measured as a function of an electrostatic force controlled by Vq. It was shown that 
Umax (Vg) can be rescaled into an universal curve that only depends on the geometrical 
parameters of the structure (L, D and H) and on Y. This scaling law, together with our 
on-chip method constitutes a very useful tool for designing actuators and in general CNT- 
based NEMS where a precise knowledge of static deflection is required. We note that the 
principle of the method used here can be extended to other geometries like that of suspended 
cantilever or multiple gate structures. 
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ported by ACI "nanosciences et nanotechnologies" programs (CNRS) and the Swiss National 
Science Foundation and its NCCR "Nanoscale Science". 
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FIG. 1: a) Schematic picture of a nanoelectromechanical doubly clamped suspended CNT three 
terminal device. The distance between the CNT and the back gate electrode is labeled as H. L 
denotes suspension length and D the external diameter of the CNT. The deviation from the straight 
line is denoted by y(x). b) SEM image of a typical device used in AFM experiments (scale bar: 
500 nm). The CNT is connected with two metallic electrodes, used both as conducting electrodes 
and anchor pads. The heavily n doped Si substrate (10 19 cm -3 ) act as a back-gate. 
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FIG. 2: Feedback voltage applied to the piezoelectric element V p i ezo as a function of the gate 
voltage Vg at two different positions on device #1 (L = 600 nm, D = 10 nm ): 1 at the 
middle of the CNT ; 2 on the sacrificial oxide of the silicon substrate (see insert). In the first case, 
increasing \ Vq\ produces the deformation of the CNT, the cantilever tip has more room to oscillate 
and the oscillation amplitude increases. To maintain constant the oscillation amplitude of the tip 
(tapping mode), a positive V p i ezo is applied by the feedback loop. Conversely, in the second case, 
increasing \Vg\ increases the attractive electrostatic force between the tip and the substrate, the 
oscillation amplitude decreases and a negative V p i ezo is applied by the feedback loop. The actual 
CNT deflection is obtained by subtracting from the curve measured at position 1 the one measured 
at position 2 previously divided by the relative permittivity of the silicon dioxide (dotted line). 



10 




0.01 



0.01 



□ 


# 1 


o 


#2 


o 


#3 


A 


#4 



0.06 



V JV) L(ijm) 2 / D(nm) 



3-/3/2 



FIG. 3: a) Maximum deflection Umax as a function of Vq for devices #1 (L = 600 nm, D = 10 
nm ), #2 (L = 770 nm, D = 13.5 nm ) and #3 (L = 480 nm, D = 10 nm ). Dash line: guide 
to the eye that has Vq dependence, b) Umax as a function of Vq rescaled using Eq|3| (see text) 
for four different devices: #1,#2,#3 and #4 {L = 730 nm, D = 15 nm ). Insert: Calculated 
Young's modulus obtained from the scaling constant K (Eq0J for the four devices. Error bars are 
estimated from the uncertainty on L and D. 
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FIG. 4: a) Calculated maximum deflection Umax as a function of Vq for different L, D parameters: 
a (L = 1200 nm, D = 13 nm), b (L = 1000 nm, D = 13 ran), c (L = 800 nm, L» = 10 ran), 
d (L = 800 nm, D = 15 nm), e (L = 800 nm, D = 20 nm), f (L = 500 nm, £> = 13 nm), 
g (L = 800 nm, D = 8 nm). b) Same results rescaled according to Eq|l] (see text). Dash line: 
Analytical result obtained for induced stress T = 0. 
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